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Abstract 

We give a complete characterization of the possible response matrices at a fixed 
frequency of n-terminal electrical networks of inductors, capacitors, resistors and 
grounds, and of n-terminal discrete linear elastodynamic networks of springs and 
point masses, both in the three-dimensional case and in the two-dimensional case. 
Specifically we construct networks which realize any response matrix which is com- 
patible with the known symmetry properties and thermodynamic constraints of 
response matrices. Due to a mathematical equivalence we also obtain a character- 
ization of the response matrices of discrete acoustic networks. 
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1 Introduction 

It is well known that composites built from high contrast constituents can have mod- 
uli or combinations of moduli which are not usually seen in nature. For example, by 
combining stiff and compliant phases one can obtain composites with a negative Pois- 
son's ratio, having a high shear modulus but low bulk modulus [U [2]. More generally 
one can construct anisotropic composites having any desired positive definite elasticity 
tensor 0, HI IS]. Composites have recently been constructed with a negative refractive 
index, having a negative electrical permittivity and a negative magnetic permeability 
over some frequency range [6] . They have also been constructed with a negative effective 
density and with a negative effective stiffness 0, IE] over a range of frequencies. Less well 
known, though perhaps more interesting, is the fact that the equations describing the 
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macroscopic behavior of composites built from high contrast constituents can be entirely 
different from those seen in nature. For example one can obtain materials with macro- 
scopic non-Ohmic, possibly non-local, conducting behavior, even though they conform to 
Ohm's law at the microscale [9j [TOj [HI [T3], [TJ], materials with a macroscopic higher 
order gradient or non-local elastic response even though they are governed by usual linear 
elasticity equations at the microscale [T5l [T6l 15]). materials with non-Maxwellian macro- 
scopic electromagnetic behavior [17], even though they conform to Maxwell's equations 
at the microscale, and materials with macroscopic behavior outside that of continuum 
elastodynamics even though they are governed by continuum elastodynamics at the mi- 
croscale [18J. It is becoming increasingly apparent that the usual continuum equations 
of physics do not apply to materials with exotic microstructures. 

One would really like to be able to characterize the possible macroscopic continuum 
equations that govern the behavior of materials, including materials with exotic mi- 
crostructures. A strategy for doing this was developed by Camar-Eddine and Seppecher 
[131 [5]). Basically, the idea is to first show that one can use a continuum construction to 
model a discrete network, consisting of nodes (terminals) which are strongly coupled to 
the continuum matrix and other nodes that are effectively hidden because they occupy 
vanishingly small volume and are essentially uncoupled with the continuum matrix [alter- 
natively, following the ideas of Milton and Willis [19] these nodes might be in an region 
of the material that is declared to be hidden, where the behavior of the fields do not 
influence the chosen macroscopic descriptors]. The next step is to characterize the possi- 
ble responses of discrete networks, in which one only is interested in the behavior at the 
terminals. The final step is to characterize the possible continuum limits of these discrete 
structures. This program was successfully carried out for three-dimensional conductivity 
[T3] and three-dimensional linear elasticity [5], giving a complete characterization of the 
possible macroscopic equations, under some assumptions such as that the source term 
does not vary on the microscale, and that the macroscopic descriptor is a single potential 
(for electrical conductivity) or a single displacement field (for linear elasticity). 

Our ultimate goal would be to characterize the possible macroscopic electrodynamic, 
acoustic, and elastodynamic equations, achievable under the assumption that the mi- 
crostructure does not vary with time, and also when this assumption is relaxed. A 
more reachable objective would be to characterize the macroscopic behavior under the 
assumption that the fields are time harmonic, oscillating at a fixed real frequency u. 
This paper is devoted to such a characterization, for discrete dynamical electric net- 
works with grounds, discrete acoustic networks, and discrete elastodynamic networks, 
anticipating that this will be key to understanding the possible macroscopic limits in 
continuum systems. Curiously, the characterization of the response tensors in the dy- 
namic case turns out to be easier than in the static case. In the static case, the possible 
response tensors of n-terminal resistor networks in three-dimensions was essentially char- 
acterized by Kirchhoff and is known as the generalized Y — A theorem: any ra-terminal 
network is equivalent to an n-terminal network having no internal nodes and with up to 
n(n — l)/2 resistors connecting the terminal pairs. However, to our knowledge, there is 
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no such characterization in two- dimensions. One exception is for circular planar resistor 
networks, where the terminals are at the boundary of a circle, and the network is con- 
tained within the circle. For this class of planar network Curtis, Ingerman, and Morrow 
[20] have completely characterized the possible response matrices. For static n-terminal 
spring networks Camar-Eddine and Seppecher [5] obtained a complete characterization 
of the possible response matrices in three-dimensions, but again the two-dimensional case 
remains an open problem. 



2 Electrical Circuits 

2.1 The lossless electrical case 

To begin with, let us treat the case of an n-terminal network consisting only of capacitors 
and inductors. An n-terminal network is a set of n + m nodes P r . Each pair (P r ,P s ) 
of nodes may be connected by capacitors and/or inductors. The n first nodes, called 
the terminals of the network, are connected to the exterior. When the terminals Pi, 
. . . ,P n are respectively submitted to voltages V\e~ wt . . . ,V n e~ lu/t , the complex currents^ 
entering the n terminals take the form iAie~ lult / 'u , . . . , iA n e~ luJt / 'u . If we denote, in the 
same way, iI rtS e~ tuJt /uj the complex current flowing to node r from node s, the linear 
behavior of the capacitors and inductors connecting the two nodes leads to the relation 

/ rs Is,r k r}S (V s V^.) (2.1) 

where the coefficient k r s is k r s — 1/L for a single inductor while fc r , s = —w 2 C for a single 
capacitor, where L is the inductance and C is the capacitance. Of course, fc r , s = when 
the two nodes are not connected. So any constant k r>s G IR is possible. Note that, in this 
description different nodes simply joined by wires are considered as a single node. We 
do not allow the terminals to be in this situation (no short circuit). 

When an internal node P r (with r > n) is not connected to the ground, Kirchhoff 's 
current law must apply. We have 

n+m 

^ = °> ( 2 - 2 ) 

s=l 

(in which we set I r ^ r = 0). At each terminal P r (with r <n) the same law reads 

n+m 

A + ^ = °- (2.3) 

s=l 

But an internal node can be connected to the ground. At such a grounded node P r , 
a current can flow toward the ground. Equation (12.21) does not apply anymore and has 

1 We have chosen to keep track of the parameters A r rather than the currents to unify the mathematics, 
and make the connections with the discrete elastic models discussed in this paper more transparent. 
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to be replaced by 



V r = 0. 



(2.4) 



We only consider circuits for which u is not a resonance frequency. Then the response 
A = (v4 r )" =1 depends in a linear way on the applied voltages V = (V^.)" =1 : there exists 
an n x n matrix W with real coefficients W r<s such that 

n 

A r = J2w r , s V s . (2.5) 

s=l 

It is well known that this matrix is symmetric 

W r , s = W s>r , (2.6) 

being the Schur complement of the matrix characterizing the response when all nodes 
are regarded as terminals, which is clearly symmetric. 

Our goal is to characterize the set of matrices W which can be obtained as a response 
matrix of a general (grounded) network but we will also consider two possible restrictions 
for the networks : 

Ungrounded networks : In this case, the network is not connected to the ground and, 
at each internal node P r (with r > n), equation (12.21) applies. In that case the response 
to a uniform voltage (Vi = V2 = . . . = V n ) is zero and the matrix W has to satisfy 

n 

Vr, £V P ,. = 0. (2.7) 

s=l 



Special grounded networks : For reasons which will become clear in section [2~4| where 
we treat acoustic networks, we pay particular attention to circuits in which inductors are 
used only to join ungrounded nodes while capacitors are only used to join an ungrounded 
node to a grounded one. Owing to (12. 4p the grounded nodes are easily eliminated in a first 
step when computing the response matrix of the circuit and at any node P r connected to 
a grounded node with a capacitor with capacitance C r , equation (12.21) has to be replaced 
by 

n+m 

h,s = -uj 2 C r V r . (2.8) 

8=1 

We say that P r "has capacitance C r " . Such circuits, we call "special grounded networks" , 
can then be considered as networks of nodes P r with capacitance C r only joined by 
inductors. 

The cases of grounded and ungrounded networks are very similar. Indeed, when con- 
sidering an ungrounded network, we can assume, without loss of generality, that one of 
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Figure 1 : Example [T] 

the terminals, let say P n , has voltage and decide to call it the "ground" . Due to the con- 
straint (12.71) the response matrix of the network will be determined by the (n — 1) x (n— 1) 
reduced response matrix W = (W rtS )r^=% where one deletes the n-th row and column 
from W. Considering P n as an internal node instead of a terminal, transforms the n- 
terminal ungrounded network in a (n— l)-terminal grounded one. The response matrix of 
this network coincides with the reduced matrix of the initial network. Reciprocally, when 
considering a grounded network, it suffices to connect all the grounded nodes together, 
making so a single node, and to consider this node as a new terminal. We then obtain 
an ungrounded new network, the reduced matrix of which corresponds to the response 
matrix of the initial network. The problems of finding all possible response matrices 
for grounded or ungrounded networks are identical as far as there are no topological or 
physical restrictions preventing from connecting together all the grounded nodes. 

Let us now consider some simple examples of special grounded n-terminal networks. 
Of course, the same response matrices may be obtained more directly as the response 
of general grounded networks. We also leave the reader to construct the ungrounded 
networks with the same reduced response matrix. 

Example 1 Let k e H* (the set of non-zero reals) and consider the simple one-terminal 
network in which the terminal is connected to an internal node of capacitance C = 
pj|p > with an inductor of inductance L — ||r (see fig. [!]). The response ma- 
trix W is that of the inductor in series with the capacitance: 

W = ([L-l/(^C)}-') = (k). (2.9) 

Using copies of this circuit with k > in a network is equivalent to allowing the 
use of internal nodes with negative "capacitance" when constructing special grounded 
networks. 

Example 2 Let k e R* and set C : = 2ig±* ; L := jL, C := ^+^. We consider the 
following two-terminal network (n = 2): terminals 1 and 2 have capacitance C. They 
are joined to an internal node P 3 with two inductors of the same inductance L. P 3 has 
capacitance C (see fig. [||). We let the reader check that the response matrix of this circuit 
is the very elementary matrix 




(2.10) 
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Figure 2: Example [2] and [3] 



Example 3 If we modify only the value of C in the previous example by setting C := 
2 ^ k \+ k "> ^ we g e t foe response matrix 



When k < this response can more directly be obtained by a unique inductor with induc- 
tance —k' 1 joining the two terminals. When k > the circuit is equivalent to a unique 
capacitor with capacitance ku~ 2 joining the two terminals. 

So, starting with a general grounded network and replacing all capacitors by copies 
of this circuit leads to a special grounded network with the same response matrix. The 
restriction to special grounded networks does not reduce the set of possible response 
matrices. 

Superposition principle : We assume that all the components (capacitors, inductors, 
wires and nodes) the network is made of occupy arbitrarily small volume. Then, in the 
three-dimensional case, the physical placement of the components has no importance and 
no crossing problem occurs. When considering two networks sharing the same terminals 
(Pi,..., P n ), by using (if necessary) a suitable distortion of one of the networks, we 
can assume that the internal components of the two networks do not intersect. So 
the response matrix of the network obtained by superposition is simply the sum of the 
response matrices of the two initial networks (because the two networks share a common 
set of voltages at the terminals, and the current flowing into each terminal is a sum of 
the currents flowing into each of the two networks through that terminal). Note that 
any n-terminal network can be considered as a m-terminal network (with m > n) in 
which m — n terminals are not connected. This superposition property shows also that, 
in dimension three, the problems of finding all possible response matrices for grounded 
or ungrounded networks are equivalent. 
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Then we easily get the following 



Theorem 1 In three- dimensions at a fixed given frequency, any real symmetric matrix 
S can be realized as the response matrix of a special grounded network. It can also be 
realized as the reduced response matrix of an ungrounded network. 

Proof: Let n be the dimension of the matrix. Let us construct a special grounded 
network the response matrix of which is S. We first consider the superposition of n 
copies of Example [U The constant k used in the copy attached to terminal r is chosen by 
setting k = S r ^ r . So the response matrix of the superposition coincides with the diagonal 
part of S. Then we superimpose on the previous network n( -"~ 1 ^ copies of Example [2J 
The constant k used in the copy attached to the pair of terminals (P r , P s ) is chosen by 
setting k = SV, S . This fixes the off-diagonal elements of the response matrix. 

An ungrounded network, the reduced response matrix of which is S, can be obtained 
using the correspondence we already described. □ 

Note that for ungrounded networks a much simpler construction is possible. Given a 
real (n + l)-dimensional symmetric matrix S, whose row sums are zero, we can realize S 
as the response matrix of an ungrounded (n + l)-terminal network by connecting every 
pair of terminals (P r ,P s ) with a component with constant k = —S r>s : then all the off- 
diagonal elements take their desired values and the diagonal elements automatically take 
the correct values by the constraint (12. 7p . 

2.2 The lossy electrical case 

Now let us extend our definition of networks by allowing resistors in the connections 
between nodes. The only change in our analysis is the fact that the constant k in 
equation (12.11) is no longer real. Indeed for a single resistor with resistance R connecting 
nodes r and s the relation (12. ip holds with k = —iuj/R. More generally k has a negative 
imaginary part. 

We also slightly extend the definition of special grounded networks by allowing any 
ungrounded node to be joined to a grounded one by a resistive capacitor. The induc- 
tors could also be resistive, but in our constructions we will still require that pairs of 
ungrounded nodes be connected only by perfect inductors. 

The response matrix W of such circuits is complex, symmetric, with negative semidef- 
inite imaginary part, 



This well-known constraint reflects the second law of thermodynamics that the circuit 
can transform electrical energy into heat but not the reverse. To see this directly it is 
easy to check that (12.121) is satisfied for a circuit in which all nodes are terminals, and as 
a result the quantity 



(ImV) • (ReA) - (ReV) • (ImA) = -(ReV) ■ ImW(ReV) - (ImV) ■ ImW(ImV) (2.13) 



ImW < 0. 
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(which is proportional to the time averaged power dissipation) is always non-negative, 
where V = (Vi, V 2 , . . . , V n ) and A = (Ai, A 2 , . . . , A n ) = WV. This remains true if some 
of the A r are zero which corresponds to a network with internal nodes. Then the left hand 
side of (12.13p is just a sum involving the V r and A r at the terminals and the algebraic 
identity implies (12.121) holds for the response matrix W of a network with internal nodes. 
It is easy to check that the response matrix has the property ( 12. 12ft if and only if the 
reduced response matrix has a negative semi definite imaginary part. 
We have 

Theorem 2 In dimension three, every symmetric complex matrix S with negative semi- 
definite imaginary part is realizable as the response matrix of some special grounded 
network. It is also realizable as the reduced response matrix of an ungrounded network. 

Proof: Let us first reconsider Example [U Let k be a complex with negative imaginary 
part and fix L — t|j (which is still a positive real and corresponds to a perfect inductor) 

and C = ( 2 fc-|fc|)o; a (which has positive real and imaginary parts and then corresponds to 
a resistive capacitor). The response matrix is again W = ( k). Thus we are allowed to 
use internal nodes with any complex "capacitance" with positive imaginary part when 
constructing special grounded networks. 

Now, let us consider the real and imaginary parts of S = S re + iS zm . They are n x n 
symmetric matrices with real coefficients. Thus we introduce the n eigenvalues (k m )^ n=l 
of S™. As S* m is a negative semidefinite symmetric matrix these eigenvalues (k m ) are 
non-positive reals, and the associated eigenvectors can be chosen to be orthonormal. We 
denote by ( ) the n-component eigenvector associated with the eigenvalue 

k m . 

Owing to Theorem [T] we know that there exists a lossless electrical circuit with 2n 
terminals with the real response matrix S with entries S r>s defined by 

S r>s = S'rs) if r <n and s < n, S r>s = 0, if r > n and s > n, (2-14) 

S r ,n+m = Sn+m,r — ~k m a™, if t < 7i and m G {1, . . . n}. (2-15) 

Let us endow each terminal P n + m (for m G {1, . . . n}) with the effective capacity C n+m = 
—ik m uj~ 2 and consider all these terminals as internal nodes. At each node P n + m (for 
m G {1, • • • n}) we have 

n 

k m J2<V S = -uj 2 C n+m V n+m = %k m V n+m (2.16) 

s=l 

and, at each terminal P r , for r < n: 

n n 

A r = J2 SZV, + ~ km < V n+m- (2-17) 
s=l m=l 
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Using the first equation to eliminate the terms involving V n+m in the second equation, 
we get 



2.3 The construction in two-dimensions 

If we think of inductors as coils of wires, then it does not make much physical sense 
to consider a planar circuit. However metals such as gold or silver or other plasmonic 
materials can have an electrical permittivity which is close to being real and negative 
over certain frequency ranges and a rectangular block of such a material can function as 
an inductor [2TJ 122] . 

Now the cases of grounded and ungrounded circuits are quite different. In the first 
case the ground is freely distributed at any internal node while in the second case only 
the nodes which can be connected with a fixed terminal may be grounded. In two- 
dimensions we have the important topological restriction that no two edges are allowed to 
cross without intersecting at a common node. This would suggest that the superposition 
principle does not apply. Surprisingly, we still have the following 

Theorem 3 Every symmetric complex matrix S with negative semidefinite imaginary 
part is realizable as the response matrix of some planar special grounded network. It is 
also realizable as the reduced matrix of a planar ungrounded network. 

Note that the circuits described in Examples [Q, [2] or [3] are planar circuits. Let us add 
a new example. 

Example 4 We consider a planar four-terminal ungrounded network. The four termi- 
nals are numbered clockwise 1, 2, 3 and 4- Let k e IR*. We join the pairs of terminals 
(Pi^Py), (P2,P$), (P3, P4) and (Pi,Pi) by four identical components with constant —k. 
We introduce an internal node P5 and join each terminal to P5 by components with con- 
stant 4k. At any terminal r 6 {1,2,3,4}, denoting r the opposite terminal, equation 





Thus we get the desired response matrix. 



□ 



( TOj) reads 



1 





s=l 



and at node P? 



5, equation A 2. 2)) reads 



1 




(2.20) 



8=1 
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Figure 3: Example HJ the virtual crossing when k > 0. When k < capacitors and 
inductors have to be exchanged. 



From these two equations we deduce 



A r = k ^2(V a - V r ) -kJ2( V s- K) = k{Vf - V r 



(2.21) 



s=l 



The response matrix W is 



W = k 



( 


-1 





i 


o \ 







-l 





1 




i 





-1 





\ 





i 





-1/ 



(2.22) 



Note that, when k is very large this four-terminal circuit is an approximation of a "virtual 
crossing". The network is equivalent to two connections with constant k joining terminals 
1 and 3 and terminals 2 and 4 without intersecting. 

The same matrix can be obtained as the response matrix of a planar four-terminal 
special grounded network (under our assumption that the grounds are allowed to be dis- 
connected from each other). Indeed it suffices to replace any capacitor by an ad-hoc copy 
of Example [3J 



Proof of theorem We only consider the case of ungrounded networks. The other case 
can be treated in a similar way. Theorem [2] provides a three dimensional ungrounded 
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Figure 4: The removal of a crossing. 



network which has the desired response matrix. A suitable distorsion transforms this 
circuit in a planar one. But the resulting network is not a true planar circuit in the sense 
that the connections between different pairs of nodes (P r , P s ), (P t , P u ) cross without any 
physical interactions. The proof will be completed by proving that any circuit with p 
crossings is equivalent to another circuit with p — 1 crossings. So a simple induction 
argument gives us a planar network without any crossing, that is a true planar circuit. 

To remove a crossing point, we use a copy of the network described in Example HI Let 
us isolate a particular crossing of two connections (P r , P s ), (Pt, P u ), the constants of which 
are denoted respectively k r<s and kt, u - We add four internal nodes P ro , P so , P to , P Uo and 
replace the two connections by four connections (P r ,P ro ), (P s ,P So ), (P t ,Pt ), (P u ,P Uo ) 
and a copy of Example 0] (with an arbitrary constant k e H* satisfying k ^ k T ^ s and 
k 7^ kt tU ) as shown in figure HI 

We choose the constants of the connections (P r ,P ro ), (P s ,P So ), (P t ,P t() ), (P U ,P U0 ) 
respectively equal to 

kr,ro 2 ( ~ ~ J , ^s,sq k rro (2.23) 

kt >t° = 2 (l tu -l) ' k ^o = Kt (2.24) 

Note that, as k e 1R*, the imaginary parts of these constants are, like the imaginary 
parts of k r ^ s and k tjU , non-positive. Checking that the response matrix is unchanged is 
straightforward. □ 
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2.4 Application to the discretized acoustic equation 

The preceeding analysis also applies directly to the discretized acoustic equation. 

A domestic water supply network is made of tubes containing a (almost) incompress- 
ible fluid and some hydraulic capacitors. Such a situation can be found also in natural 
conditions: for instance in an unsaturated porous medium. Consider a two-dimensional 
or three-dimensional network of tubes with cavities at the junctions. Each tube contains 
a segment of incompressible, non- viscous, fluid with some density, possibly varying from 
tube to tube, moving in a time harmonic oscillatory manner in response to time har- 
monic pressures at the junctions. (There could an additional time independent constant 
pressure everywhere, but this does not affect the equations). We define the entire cavity 
associated with a junction to be the cavity at the junction, plus the remaining region in 
the tubes not occupied by the incompressible fluid. Each entire cavity contains a com- 
pressible, non-viscous, massless fluid with compressibility possibly varying from junction 
to junction. The surfaces between the compressible and incompressible fluids have some 
surface energy so that the interfaces remain flat. 

When the terminals Pi, . . . ,P n are respectively submitted to pressures pie~ luJt , . . . ,p n e~ 
the complex fluid currents entering the n terminals take the form iAie~ luJt /u, 
iA n e~ lu)t / uj . We denote, in the same way, iI r}S e~ iU>t / u the complex current flowing to 
node r from node s. Let a r s be the cross-sectional area of the tube joining nodes P r 
and P s and let m r s be the mass of the fluid contained in this tube. Since the complex 
force on the fluid segment is a ns (p r — p s ) and its complex acceleration is I r>s e~ wt /a r>s , 
Newton's law of motion implies 

Ir,s = ~Is,r = kr,s\Ps ~ Pr) (2.25) 

where k r>s = a%. s /m rjS . 

Now consider an internal node P r . Due to the motions of the incompressible fluid 
segments in the tubes the volume of the associated entire cavity changes with time and 
the complex pressure p r in the cavity adjusts itself according to Hooke's law, 

n+m 

Ir,s = -C r U 2 Pr , (2.26) 

5=1 

where C r = V r /n in which V r is the volume of the entire cavity when the fluids are at 
rest, and k is the bulk modulus of the fluid in this entire cavity. When the junction is a 
terminal the sum must take into account the current entering the terminal. If we assume, 
without loss of generality, that the capacity of each terminal vanishes, we have 

n+m 

A r + T r,s = 0- (2.27) 

s=l 

Since fl2T26D and (I2T2TD are the direct analogs of equations (J2ID (£2D and 

(12.31) which describe special grounded networks, all the previous analysis and associated 



12 



Figure 5: A two-terminal discrete acoustic network. In the idealized model the four 
cavities contain compressible massless fluid, while the grey shaded fluid plugs in the 
five tubes contain incompessible fluid with some mass. The response of the network 
is measured by the movement of the two frictionless pistons, in response to the time 
harmonic forces acting on them which control the pressures in the terminal cavities. 

theorems apply. In particular, to get a negative effective "bulk modulus" in a cavity 
we just follow the approach of Fang et.al.[8] and connect that cavity to a Helmholtz 
resonator, i.e. connect it to another cavity containing compressible fluid with a tube 
containing a plug of incompressible fluid with mass chosen so that the system is above 
resonance. 

The lossy case arises when the fluid in some, or all, of the junctions has some bulk 
viscosity, so that k has some negative imaginary part, and consequently C r in equation 
(12.261) has a positive imaginary part. The lossy case also arises if the incompressible fluid 
segments have some shear viscosity so that Darcy's law implies &v jS has a complex part 
due to the fluid permeability of the tube. 

We can conclude that any response is possible for an acoustic discrete system provided 
that the total dissipation is non-negative. 

3 Discrete Elastodynamics 

We now turn our attention to networks of springs with point masses at the nodes. We 
emphasize that our analysis is for idealized linear networks and that we do not consider 
questions of stability, and in particular stability against buckling. 

Let us denote u r e~ luJt the displacement of node P r (u is a three dimensional complex 
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vector) and F T)S e~ luJt the force exerted on node P r by the spring (if any) joining P s and 
P r . In the same way, for any terminal P r (r < n), let us denote A r e~ iuJt the additional 
external force applied on the system at that terminal. At each interior node P r , Newton's 
law applies and we have 

n+m 

F r s = —m r uj 2 u r , (3.1) 

s=l 

(which is the analog of (12. 8p ) where m r denotes the mass of node P r . At a terminal P r 
(r < n), we have also to take into account the external force: 

n+m 

A r + F riS = — m r u 2 u r . (3-2) 

s=l 



3.1 The purely elastic case 

To begin with, let us treat the case of purely elastic n-terminal networks where there is 
no damping in the springs. Between each pair of nodes P r and P s , located at positions 
x r and x s that are linked by a spring, Hooke's law applies, and we have 

x — x 

F r)S = — F S)r = k r ^ s n r ^ s ® n rjS • (u s — u r ), where n rs = - — — , (3-3) 

W^-s II 

(which is the analog of (12. ip ) where k rjS = k s ^ r is the (positive real) spring constant. 

The elastodynamic response of the network is governed by a matrix W which has 
second order tensors W r s as its entries, and links the set of forces A = (A 1; A 2 , . . . , A n ) 
with the set of displacements U = (u 1; u 2 , . . . , u n ), through the relation 

n 

A r = ^W rs u s . (3.4) 

8=1 

This matrix W is real and has the symmetry property that, for any r, s in {1, . . . , n}, 

W r , s = (W s , r ) T (3.5) 

and the proof of this property is similar to the proof of (12. 6p in the electrical case. 

Example 5 The simplest non-trivial two terminal network just consists of terminals P\ 
and P2 joined by a spring with constant k and direction n = ni ;2 . According to / tff.gp the 
matrix W is 

= / fcn<8>n -£;n<g>n\ 
y — fcnCEDn kn®n J' 

Example 6 Let us consider the two terminal network in which the terminals P\ and 
P 2 are joined to three internal nodes P3, P^, P5 making two non- degenerate simplexes 
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(Pi, P3, P4, P5) and (P2, P3, P4, P5) . To each edge of these simplexes corresponds a spring. 
Nodes have no mass. In structural mechanics such a structure is called a simple truss. 
Its response vanishes when the applied displacements u 1; u 2 correspond to a rigid motion. 
So the response matrix corresponds to a non-negative quadratic form depending only on 
(ui — u 2 ) • n 12 . It takes the form Ii3. 6\) . The constant k can be tuned by multiplying all 
the constants of the truss by a common positive factor. 

Note that the choice of the position of the internal nodes is quite free. A given finite 
set of points can easily be avoided. Moreover, in dimension three, the internal nodes can 
be chosen in such a way that the five segments (Pi, P 3 ), (Pi, P 4 ) (P 2 , P3), (P 2 , P4) (P3, P4) 
do not intersect a given finite set of straight lines (but maybe at terminals Pi, P 2 ). 

Remark 1 Replacing a single spring in a network by a structure described in Example 
will not change the response matrix of the network. In dimension three, making all 
the needed replacements we can restrict our attention to networks in which any different 
springs do not intersect and have different directions. 

Example 7 Let be a non vanishing real and n a unit vector. Consider the very simple 
one-terminal spring network, where there is only one spring with constant k = k\^. = 
2^~jjjj ^ 2 Unking terminal Pi with a single interior node P 2 chosen in such a way that 

ni j2 = n. Terminal Pi has no mass while the mass of node P 2 is m = ^ . We have the 
equations 

F 2) i = k n g) n ■ (ui — u 2 ) = — muj 2 u 2 = A 1 . (3.7) 

The elimination of u 2 leads to Ai = — ^ m ^ w i n ® n ■ ui = —ficu 2 n (g) n • Ui. This system 
endows P\ with the tensor valued "effective mass" : M = /m (g> n which can either be 
positive or negative semidefinite depending on the sign of \x. The physical reason that one 
can obtain negative values of \x is that these are achieved when the spring-mass system 
is above resonance, i.e. when k — muj 2 < 0, and as a result the mass oscillates 180° out 
of phase with the motion of the terminal. 

Note that the choice of the position of the internal node is free on the straight line 
(xi, n). Thus a given finite set of points can easily be avoided. The spring can also be 
replaced using remark [T] in order to avoid intersections with a given finite set of straight 
lines. 

Now let M be any real symmetric tensor. Superimposing up to three copies of the 
previous structure choosing for \x the eigenvalues of M and for n the corresponding 
eigenvectors of M we get 

Remark 2 Any node can be endowed by any real symmetric tensor effective mass. 

Example 8 Let K be a real, n 2 be two unit vectors and xi, x 2 be two distinct points 
such that at least one of the two vectors m\\, n 2 is not in the direction x 2 — xi. We 
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consider the two terminal network consisting of terminals Pi, P 2 at points xi, X2 and 
two internal nodes P3, P4 placed in such a way that n li3 = n 1; n 2j4 = n 2 and v := n 3i4 
is neither colinear with ni nor n 2 . VFe introduce two new unit vectors Wi and w 2 which 
complete respectively the basis (n 1; v) and (n 2 ,v). 

Springs with constant k = \K\ join pairs (Pi, P3) and (P 2 , P4). ^4 spring with constant 
k' = 2\K\ — K joins (P3, P4). Nodes P3 and P4 are endowed respectively with the effective 
masses k'oj~ 2 {y (g) ni + ni ® v + wi ® wi) and /c'co>~ 2 (v ® n 2 + n 2 ® v + w 2 ® w 2 ). 

y4s noticed in Example^ such effective masses need the introduction of extra springs 
and internal nodes. Again we have a large freedom in the choice of the position of the 
nodes P 3 , P 4 on the lines (xi,ni) and (x 2 ,n 2 ) and we can avoid any given finite set of 
points. Owing to RemarkUl we can also construct this structure avoiding any intersection 
with a given finite set of straight lines. 

Let us introduce the dual basis (n*,v*,w^) o/(ni,v, wi) (i.e. satisfying n* • ni = 1, 
n^; • v = 0, • wi = 0, v* • ni = ; v* • v = 1 ; v* • wi = 0, w* • ni = 0, w] 1 • v = 
and wl ■ wi = 1) and in the same way the dual basis (n 2 , v°, w 2 ) of (n 2 ,v, w 2 ). The 
displacement of nodes P3, P4 are respectively written in the form 

u 3 = anl + ov* + cw*, u 4 = dn 2 + ev° + /w 2 . (3.8) 

At nodes P3 and P4, equation (3. 1\) reads 

— k(xii ■ ui)ni + kaxii — k'(e — 6)v = k'(a\ + 6ni + cwi) 
—k(n 2 ■ u 2 )n 2 + kdn 2 + k'(e - 6)v = k'(d\ + en 2 + fw 2 ) (3.9) 

from which we deduce b — e = a = —d = |^r(ni • Ui — n 2 • u 2 ) and c = f = 0. Now at 
terminals Pi and P 2 equation Ii3. || ) reads 

Ai = k(iii ■ Ui)ni — kaxii — c<j 2 Mi ■ ui 

A 2 = k(n 2 ■ u 2 )n 2 - A;dn 2 - co 2 M 2 ■ u 2 , (3.10) 

where Mi and M 2 denote the effective masses of terminals Pi, P 2 which we have not yet 
fixed. Then the response matrix of the network is then 

= f -cj 2 Mi + (\K\ - AT)m ® ni Kn x ® n 2 

V Kn 2 ®ni -w 2 M 2 + (|K| -iT)n 2 ®n 2 

T/ie key feature of this response matrix is that the off-diagonal matrix is proportional 
to ni ® n 2 . 27ms conW nave been anticipated since the spring joining terminals (Pi,P 3 ) 
exerts a force on Pi in the direction n 1; and this force can only depend on u 2 through the 
component of u 2 in the direction n 2 of the spring joining terminals (P 2 ,P 4 ). 

Example 9 Considering in the previous example the particular case ni = n 2 = n and 
choosing the appropriate values for the tensors Mi and M 2 we get 

f-Kn®n Kn®n\ 
VV -\ Kn®n -Kn®n) ' [6AZ) 
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which is similar to the response matrix of a single spring but where the constant K can 
be negative. More important is the fact that in this structure the direction of action n is 
no longer correlated with the direction of the vector x 2 — x x . 

Remember however that we have the restriction that n cannot be in the direction 
x 2 — Xi. But we can get rid of this restriction by considering two copies of this structure: 
one of these copies joins terminal P\ to an internal node P3 placed at a point x 3 such 
that 11^3 is not parallel to x 2 — xi while the other one joins terminal P 2 to P3. In both 
copies the constant is 2K and the direction of action is n = ni j2 . It is easy to check that 
the response matrix of such a structure is still given by $3.12\) . In that way we actually 
get a virtual spring with possibly negative spring constant. 

This makes free the position of the internal nodes: indeed, in any network, we can 
change the position x r of a node P r to any other position x! r , replacing all the springs 
joining P r to other nodes P s by a copy of Example [9] with ni = n rjS . Clearly the response 
matrix will remain unchanged. Thus we have 

Remark 3 Any network has an equivalent network the internal nodes of which avoid a 
given finite set of points. 

Combining Remark [3] with Remark [1] enables us to assume when considering two 
different networks that they do not share any internal node and that the springs of the 
different networks do not intersect. Then we have 

Remark 4 Superposition principle : In three dimensions, if Wi and W 2 are two realiz- 
able response matrices each associated with n-terminals in the same positions P\, P 2 , . . . , P, 
then the response matrix Wi + W 2 is also realizable. The network which realizes the ma- 
trixWi J rY\?2 is just a superposition of suitable modifications of the networks which realize 
Wi and W 2 . 

Example 10 Choosing in Example the appropriate values for the tensors Mi or M 2 
we can get 



As any matrix W is the sum of rank one matrices, the superposition principle implies 
that, for any matrix W, there exists a two-terminal network the response matrix of which 



And we obtain the following 

Theorem 4 In three dimensions given any set of n points xi,x 2 ,...x n , and any real 
n x n matrix W with second order tensor entries Wjj satisfying the symmetry proper- 
ties Ii3. 5\) . then there is purely elastic network with terminals Pi,P 2 , . . . P n at positions 
xi, x 2 , . . . x n and realizing W as its response matrix. 




(3.13) 



is 




(3.14) 
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Proof: It is enough to attach to each pair (P r , P s ) a copy of Example [TOl in which W is 
chosen to be W r s , then to endow each terminal P r with the effective mass corresponding 
to the symmetric matrix W rr . Then we conclude using the superposition principle. □ 

3.2 Elastodynamic networks with damping 

An elastodynamic network with damping is a network with point masses at the nodes 
and viscoelastic springs joining the nodes. If we allow viscous damping in the springs, 
the constant k in (I3.3P becomes complex with a non-positive imaginary part. (The real 
part of k is still non-negative, and masses are still non-negative reals.) Then the matrix 
W is complex and symmetric, with negative semi definite imaginary part, 

Imag>V<0, (3.15) 

which reflects the second law of thermodynamics that averaged over time the network 
can transform mechanical energy into heat, but not the reverse. The proof of (13. 15ft is 
similar to the electrical case. 

Let us revisit the examples we gave in the previous section. Examples [5] and [6] are 
unchanged: the constant k in the response matrix is now complex with a positive real 
part and a negative imaginary part. Example [7] is still valid : indeed for any complex 
jj with positive imaginary part the constant k defined by k = 2j^T\ uj2 nas a ne g a tive 
imaginary part and a positive real part. Then remark [2] can be generalized in 

Remark 5 Any node can be endowed by any complex symmetric tensor effective mass 
provided that its imaginary part is positive semidefinite. 

Theorem 5 In three dimensions, given any set of n points xi,X2, . . . x n , and any com- 
plex matrix W with second order tensor entries W rjS satisfying the symmetry properties 
\3. 5\) and the constraint Ii3.15\) . then there is an elastodynamic network with damping 
realizing W as its response matrix, 

Proof: Let us consider the real and imaginary parts of W = W re + iW im . They are 
n x n matrices with 3x3 real entries denoted respectively W£ e s and W*™ and can be 
identified with 3n x 3n symmetric matrices with real coefficients. Thus we introduce the 
3n eigenvalues {k m )^ =1 of W jm . As W im is a negative semidefinite symmetric matrix 
these eigenvalues (k m ) are non-positive reals, and the corresponding eigenvectors can 
be chosen to be an orthonormal set. We denote by (a™, a™, . . . a™) a 3n-component 
eigenvector (identified with an n-entry vector where the entries are 3-component vectors) 
associated with the eigenvalue k m and we introduce an extra unit vector b. 

Owing to Theorem S] we know that there exists a (no damping) elastic network with 
An terminals with the real response matrix W with entries W r . jS defined by 

W rjS = W^ e s , if r < n and s < n, W r>s = 0, if r > n and s > n, (3.16) 
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W r>n+m = W^ +mjr = -k m a™ ® b, if r < n and m 6 {1, . . . 3n}. (3.17) 

Then, owing to Remark[5l let us now endow each terminal P n + m (for m £ {1 ; • • • 3^}) 
with the effective mass tensor M n+m = —ik m u~ 2 1 and consider all these terminals as 
internal nodes. At each node P n+m (for m G {1, . . . 3n}) we have 

n 

fc m J](a™ • u r ) b = -u 2 M n+m ■ u n+m = tk m u n+m (3.18) 

r=l 

and, at each terminal P r , for r < n: 

n 3ra 

A r = J2 W Z - U s + J2 ~ km ( h ■ U n+m)B?- (3-19) 

s=l m=l 

Using the first equation to eliminate the terms involving u n+m in the second equation, 
we get 

n / 3n \ 

s=l \ m=l / 

Thus we get the desired response matrix. □ 



3.3 Planar elastodynamic networks 

As in the electrical case, in two- dimensions, we have the important topological restriction 
that no two edges are allowed to cross without intersecting at a common node. Now we 
have the additional restriction that a spring between two nodes must lie along the segment 
joining those two nodes. Despite these restrictions we have: 

Theorem 6 Theorems^ and{5\ still hold true for planar networks. 

Proof: Example [6] can be adapted to the planar case : the simplexes we used are now 
simply triangles. However Remark[T]is no longer valid. Due to the topological restrictions, 
Example M cannot be used to avoid crossings. It still can be used to change the direction 
of the springs in a network. So we can assume that any crossing point is a generic one, 
which means that only two springs are crossing at that point and the angle they make is 
non-zero. 

Let us allow, for a while, springs to intersect without interacting and let us call pseudo- 
planar such networks. In this setting, Example [7] and [H] and the superposition principle 
are still valid. Nothing is changed from the three dimensional case : we can construct a 
pseudo-planar network with any desired response matrix. Owing to the previous remark 
we can assume that all crossing points in this pseudo planar network are generic ones. 

Now let us consider two crossing springs (let us say connecting (Pi, P2) and connecting 
(P 3 , P 4 ) with constants k^ 2 and fc 3i4 ) in this network and let us replace these two springs 
by the following network: 
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Example 11 Let &i i2 and fc 3j4 be any positive reals (or complex with positive real part 
and negative imaginary part) and consider the ^-terminal network where the terminals 
Pi (i < 4=) are placed at points Xj such that the segments [x 1; x 2 ] and [x 3 ,x 4 ] have an 
intersection at a single point x 5 . The network has an internal node P5 at point x 5 . We 
assume that the nodes have no mass and that four springs join P 1; P 2 , P 3; P 4 to P5 with 
constants respectively equal to k\$ = k 2 ^ : = 2fci )2 and fc 3j 5 = fc 4j5 := 2fc 3)4 . We have 

Ai = F 1>5 = -2A;i i2 (n 1)2 ® n 1>2 ) ■ (m - u 5 ), A 2 = F 2j5 = -2k lj2 {n 1>2 <g> n 1>2 ) ■ (u 2 - u 5 ), 

(3.21) 

A 3 = F 3i5 = -2fc 3j4 (n 3i4 <g> n 3)4 ) • (u 3 - u 5 ), A 4 = F 4j5 = -2A; 3i4 (n 3i4 <g> n 3)4 ) • (u 4 - u 5 ), 

(3.22) 

F 1>5 + F 2i5 + F 3>5 + F 4>5 = 0. (3.23) 

Owing to the geometrical assumptions (n li2 ,n 3 4 ) makes a basis and we introduce its dual 
basis (n* 2 , rig 4 ). Writing u 5 = anj 2 + bn^ 4 , the previous system of equations becomes 

Ai = -2&i i2 (iii >2 • Ui - a)ni j2 , A 2 = -2k lt2 (n li2 ■ u 2 - a)ni j2 (3.24) 

A 3 = -2A; 3i4 (n 3i4 ■ u 3 - 6)n 3j4 , A 4 = -2/c 3i4 (n 3i4 ■ u 4 - 6)n 3j4 (3.25) 
nj.,2 ■ (ui + u 2 ) = 2a, n 3>4 • (u 3 + u 4 ) = 2b (3.26) 
The elimination of u 5 (i.e. of a and b) in this system leads to 

Ai = -A 2 = fci )2 (ni >2 ® n lj2 ) ■ (ui - u 2 ), A 3 = -A 4 = k 3A (n 3A ® n 3j4 ) • (u 3 - u 4 ) 

(3.27) 

The response matrix of this network is equivalent to the response of two springs joining 
directly and independently Pi to P 2 and P 3 to P 4 with constants &i j2 and A; 3j4 . 

Replacing two crossing springs by a copy of Example [TT1 removes a crossing point (and 
does not create any new one). Hence we can successively remove all crossing points in 
the pseudo planar network and obtain a true planar network with the desired response 
matrix. This analysis is valid in both purely elastic and damping cases. □ 
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